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Abstract 

In the beginning of spring, floods are the most important geomorphic hazards in Iran, destructing 
propertiesas and human lives. Aland Chai basin, located in Khoy County (northwest Iran), is also 
known as one of the basins with high potential for flood hazard due to its special geographical 
situation. This study tried to model spatial variation in flood hazard susceptibility in this basin 
using the ensemble model, FURIA-GA-LogitBoost. For this purpose, 13 effective parameters of 
flooding including lithology, soil hydrological groups, NDVI, land use, slope, aspect, elevation, 
distance to the river, river density, precipitation, topographic wetness index, stream power index, 
and sediment transport index were used. WEKA software was used to implement the research 
model and the final flood hazard susceptibility map was prepared. The study found that 
downstream areas of the basin have a high flood hazard susceptibility. These areas contain the 
most important human settlements (Khoy city) and agricultural lands and flood as a geomorphic 
hazard can seriously damage them. Considering the ROC curve and area under the curve (AUC), 
it was found that the FURIA-GA-LogitBoost model performed well in the preparation of flood 
hazard susceptibility map with coefficients of 0.861 and 0.895, respectively, in training and 
validation data. 
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جغرافیا و مخاطرات محیطی. سال دوازدهم. شمارة چهل و پنجم. بهار ۱۶۰۲ صص ۱-۲۶ 


مقاله پژوهشی 


مدل‌سازی حساسیت خطر وقوع سیل در حوضه ppl‏ الندچای بر پایه یک رویکرد طبقه‌بندی ترکیبی 
نوین (FURIA-GA-LogitBoost)‏ 


® توحید رحیم‌پور" - دکتری ژئومورفولوژی, دانشگاه تبریزه تبریزه ایران 
محمدحسین رضائی مقدم - استاد گروه ژئوموفولوژی» دانشگاه تبریز» تبریز» ایران 
سید اسداله حجازی - دانشیار گروه ژئومورفولوژی» دانشگاه تبریز» تبریز ایران 


خلیل ولیزاده کامران - استاد گروه سنجش‌ازدور و GIS‏ دانشگاه تبریز» تبریزه ایران 


تاریخ دریافت:  ۱:۰۰/۹/۲۲‏ تاریخ بازنگری: ۱۶۰۰/۱۱/۳۰ تاریخ تصویب: ۱۶۱۰/۱۲/۸ 


حکیده 


با شروع فصل بهار سیلاب‌ها به عنوان مهم‌ترین مخاطره ژئومورفیک در سطح کشور مطرح می‌شوند که 
خسارت‌های جانی و مالی فراوانی را به بار می‌آورند. حوضه آبریز الندچای واقع در شهرستان خحوی و 
شمال غرب کشور نیز به دلیل موقعیت خاص جغرافیایی جزو حوضه‌های با پتانسیل بالای خطر وقوع 
سیل شناخته می‌شود. هدف از پژوهش حاضر مدل‌سازی تغییرات فضایی حساسیت خطر وقوع سیل در 
این حوضه با استفاده از مدل ترکیبی نوين FURIA-GA-LogitBoost‏ می‌باشد. به همین منظور از ۱۳ 
پارامتر مژثر در وقوع سیل شامل لیتولوژی, گروه‌های هیدرولوژیکی خاک شاخص پوشش گیاهی, کاربری 
اراضی, ارتفا شیب جهت شیب. فاصله از آبراهه. تراکم آبراهه بارش شاخص رطوبت توپوگرافیک» 
شاخحص قدرت آبراهه و شاخص حمل رسوب استفاده شده است. جهت اجرای مدل تحقیق از نرم‌افزار 
WEKA‏ استفاده شده و نقشه نهایی حساسیت خطر وقوع سیل تهیه گردید. یافته‌های پژوهش نشان 


۱ نویسنده مسئول: ۱۹۳۹۲۸۶۹۹۸۲ Email: Rahimpour1990 ۵ gmail.com‏ 
نحوه ارجاع به این مقاله: 

رحیم پور» توحید؛ رضائی مقدم. محمدحسین؛ حجازی» سید اسداله؛ ولیزاده کامران» خلیل. ۱۶۰۲. مدل‌سازی حساسیت خطر 
وقوع سیل در حوضه آبریز الندچای بر پایه یک رویکرد طبقه‌بندی ترکیبی نوین (FURIA-GA-LogitBoost)‏ جغرافیا و 


محاطرات محیطی. ۱(۱۲). صص ۱-۲۶ 
https://doi.org/10.22067/geoeh.2022.74170.1141‏ 


۲ جغرافیا و مخاطرات محیطی شمارة اول 


می‌دهد مناطق پایین‌دست حوضه حساسیت بالایی را از نظر حطر وقوع سیل دارند. این مناطق محل تمرکز 
مهم‌ترین اجتماعات انسانی حوضه آبریز (شهر خوی) و زمین‌های کشاورزی و باغات است که سیلاب 
به‌عنوان یک مخاطره ژئومورفیک» تهدید جدی برای این مناطق محسوب می‌شود. بررسی میزان دقت نقشه 
نهایی با استفاده از منحنی ROC‏ و سطح زیر منحنی OLS (AUC)‏ داد که مدل به کار رفته در تحقیق به 
ترتیب با ضرایب ۰/۸۰۱ و ۰/۸۹۵ از نظر داده‌های آموزشی و اعتبارسنجی از عملکرد خوبی در تهیه نقشه 
حساسیت خطر وقوع سیل برخوردار بوده است. 

کلیدواژه‌ها: سیلاب. FURIA‏ الگوریتم ژنتیک, مدل LogitBoost‏ حوضه آبریز الندچای. 


۱- مقدمه 


ارائه تعریفی از سیل با توجه به ماهیت پیچیده آن. کار دشواری است و به همین دلیل تعاریف مختلفی از این 
پدیده ارائه شده است. چارلتون ' در کناب خود تحت عنوان مبانی ژومورفولوژی رودنحانه‌ای تعریف خود از سیل را 
به این صورت بیان می‌کند؛ سیل عبارت است از جریان نسبتاً زیادی که بیش از ظرفیت JOE‏ رودخانه است. 
درحالی که جریان‌های Wale‏ درون کانال اتفاق می‌افتند. جریان‌های تند دوره‌ای از سواحل کانال جریان می‌گذرند و 
روی دشت‌های سیلابی اطراف جریان پیدا می‌کنند (ثقفی و رضائی‌قدم. IVAN‏ در سرتاسر جهان, بیش از یکسوم 
سطح زمین در معرض سیلاب‌ها قرار دارده جایی که بیش از ۷۰ درصد از جمعیت در آن سکونت دارند (آکسوی و 
همکاران » ۲۰۱۹). به‌طورکلی سیلاب‌ها بر اساس وقوعشان به ٤‏ نوع طبقه‌بندی می‌شوند که عبارتند از: سیلاب‌های 
ناگهانی» سیلاب‌های رودخانه‌ای. سیلاب‌های ساحلی و سیلاب‌های شهری OS)‏ بوی و همکاران" ۲۰۱۹؛ کاستاج و 
همکاران " ۲۰۲۰). در این طبقه‌بندی سیل‌های ناگهانی از مخرب‌ترین نوع سیلاب‌ها محسوب می‌شوند که خسارات 
جانی و مالی زیادی را به دنبال دارند. سیلاب‌ها می‌توانند تأثیر عمیقی بر روی اکوسیستم‌ها و زندگی انسان‌ها بگذارند 
(الکساندر* و همکاران» ۲۰۱۱)؛ به‌طوری‌که بین سال‌های ۱۹۹۵ و ۲۰۱۵ حدود ۱۰۹ میلیون نفر از سیلاب‌ها متأثر 
شده‌اند (آلفیری" و همکاران. ۲۰۱۷). روش‌های کمی و کیفی متعددی در ادییات علمی برای مدل‌سازی سیلاب در 
Jeb‏ حوضه آبریز وجود دارد (محمود و رحمان, ۲۰۱۹). توفیق‌لو اسلام" و همکاران (۲۰۲۱) مدل‌های به کار رفته 
جهت ang‏ نقشه‌های حساسیت سیل‌خیزی را با بررسی تحقیقات مرتبط با سیل به ٤‏ دسته طبقه‌بندی کرده‌اند که 


شامل ۱- مدل‌های مبتنی بر دانش کارشناس مثل مدل فرآیند تحلیل سلسله مراتبی (AHP)‏ ۲- مدل‌های مبتنی بر آمار 


1 Charlton 

2 Aksoy et al 

3 Tien Bui et al 
4 Costache et al 
5 Alexander 

6 Alfieri 

7 Towfiqul Islam 


از قبیل نسبت فراوانی CFR)‏ ارزش اطلاعات ضریب اطمینان. رگرسیون لجستیک. وزن شواهد. منطق فازی» منطق 
فازی عصبی ۳- الگوریتم‌های یاد گیری ماشین و ۶- مدل‌های هیدرولوژیکی از قبیل SWAT‏ و HEC-RAS‏ هستند. 

تحقیقات خوبی در زمینه مدل‌سازی حساسیت خطر وقوع سیل در داخحل و خارج از کشور انجام گرفته است که 
از جمله تازه‌ترین آن‌ها می‌توان به موارد زیر اشاره کرد: روستایی و همکاران (VAT)‏ در تحقیقی به تهیه نقشه 
پهنه‌بندی سیلاب حوضه آبریز نکارود با استفاده از مدل SCS-CN‏ و GIS/RS‏ پرداختند. نتایج تحقیق OLE‏ داد که 
نواحی با پوشش گیاهی کم همراه با شیب VE‏ و کاربری اراضی کشاورزی در دامنه‌های حوضه آبریز: جزو مناطق با 
سیل‌خیزی خیلی زیاد و زیاد را در بر می‌گیرد. میرموسوی و اسمعیلی (VES)‏ در تحقیقی با استفاده از سامانه 
اطلاعات جغرافیایی و سنجش‌ازدور نواحی سیل‌خیز شهرستان داراب را پهنه‌بندی کردند. در این تحقیق از ۱۲ پارامتر 
مؤثر در سیل‌خیزی شامل شیب جهت شیب بارش ارتفاع رواناب» وضعیت شماره منحنی MCN)‏ جهعت جریان 
تجمعی» وضعیت گروه‌های هیدرلوژیکی خاک کاربری اراضی. تراکم پوشش گیاهی؛ تراکم شبکه زهکشی» وضعیت 
ارتفاعی حوضه و فاصله از آبراهه‌ها استفاده شد. نتایج تحقیق نشان داد که بر اساس نقشه‌های نهایی حطر وقوع سیل 
با دوره‌های با زگشت 0 ۱۵ ۲۵ و ۵۰ ساله» به ترتیب ٩/۲٤‏ ۱۰/۶ ۱۱/۹۰ و ۸7۵ درصد از مساحت شهرستان 
دارای خطر سیل‌خیزی خیلی زباد می‌باشد. پانت" و همکاران (۲۰۲۰) با استفاده از پارامترهای مورفومتریک و 
مورفوتکتونیک نقشه پهنه‌بندی حطر وقوع سیل و فرسایش خاک در حوضه آبریز آلاک ناندا واقع در کشور 
هندوستان را angi‏ نمودند. در این تحقیق حوضه آبریز به ۰ زیرحوضه تقسیم شده و تعداد ۵۰ پارامتر مورفومتریک و 
7 پارامتر مورفوتکتونیک مورد بررسی قرار گرفت. نتایج تحلیل آماری پارامترهای مورفومتریک و مورفوتکتونیک 
نشان داد که منطقه از نظر تکتونیک فعال بوده و بسیار مستعد وقوع سیل و متعاقب آن فرسایش خاک می‌باشد. 
توفیق‌لو اسلام و همکاران (VHT)‏ در پژوهشی با استفاده از مدل‌های پیشرفته یادگیری ماشین اقدام به مدلسازی 
حساسیت سیل در زیرحوضه تیستا واقع در شمال بنگلادش نمودند. در این تحقیق از دو مدل ترکیبی جدید 
Dagging‏ و Random Subspace‏ همراه با مدل‌های شبکه عصبی مصنوعی. جنگل تصادفی و ماشین بردار پشتیبانی 
استفاده شد. نتایج تحقیق نشان داد که مدل Dagging‏ نتایج بهتری نسبت به plu‏ مدل‌ها داشته است. 

بنابراین در تحقیق حاضر تلاش شده است با استفاده از مدل ترکبی نوین FURIA-GA-LogitBoost‏ نقشه 


حساسیت خطر وقوع سیلاب در سطح حوضه آبریز الندچای واقع در استان آذربایجان غربی تهیه گردد. 


منطقه مورد مطالعه 
منطقه مورد مطالعه این تحقیق حوضه آبریز الندچای می‌باشد که ازنظر موقعیت سیاسی در استان آذربایجان غربی 


و در شمال غرب ایران واقع شده است. این حوضه با مساحت ۰ کبلومترمربع از نظر موقعیت جغرافیایی بین 


1 Pant 
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مختصات ۱۶۳ ۳۰ ۳۸ تا ۲۲۳۳ EN‏ ۳۸۳ عرض شمالی و ۱۳ ۱۵ EE‏ تا 107۰۱۰۳ طول شرقی واقع شده 
است (شکل ۱). تغییرات ارتفاعی حوضه از ۱۰۹۳ متر در بستر رودخانه الندچای در خروجی حوضه تا ۳۱۳۸ متر 
ارتفاع از سطح دریا در کوه اورین می‌باشد. این حوضه یکی از زیر حوضه‌های حوضه آبریز رود ارس محسوب 
می‌شود که آب‌های سطحی آن پس از پیوستن به رودخانه بزرگ قطورچای به رود ارس می‌ریزد (رضائی‌مقدم و 
همکاران ۱۳۹۹). حوضه آبریز الندچای ساختاری کوهستانی دارد» به‌طوری‌که بیش از We‏ درصد مساحت آن را 
نواحی کوهستانی و مرتفع تشکیل می‌دهد. کوه اورین با ارتفاعی نزدیک به ۳۹۳۸ متر مهم‌ترین عارضه توپوگرافی در 
claw‏ ناهمواری‌های دانعل حوضه آبریز الندچای می‌باشد. آب‌وهوای این حوضه phy‏ تحت SU‏ توده‌های هوای 
غربی قرار دارد که منشاً رطوبت و بارش در فصول پایین زمستان و بهار است. 
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شکل ۱ - موقعیت جغرافیایی منطقه مورد مطالعه در استان آذربایجان غربی‎ 
مواد و روش‌ها‎ -۲ 
cal به‌منظور مدل‌سازی و تهیه نقشه حساسیت خطر وقوع سیل در حوضه آبریز الندچای گام‌های زیر طی شده‎ 


که در ادامه به توضیح هر یک از مراحل پرداخته می‌شود: 


گام اول: تعیین پارامترهای مؤثر در وقوع سیل 

در گام نخست ابتدا با بررسی پژوهش‌های انجام گرفته در رابطه با موضوع تحقیق و همچنین داده‌های موجود در 
منطقه موردهطالعه, ۱۳ پارامتر مؤثر در حساسیت خطر وقوع سیل انتخاب شد. پارامترهای مورداستفاده در این تحقیق 
شامل لیتولوژی, گروه‌های هیدرولوژیکی SE‏ شاخص پوشش گیاهی» کاربری اراضیء شیب جهت شیب ارتفاع» 
فاصله از آبراهه تراکم آبراهه بارش شاخص رطوبت توپوگرافیکه شاخص قدرت آبراهه و شاخص حمل رسوب 
ela‏ 

گام دوم: برآورد همبستگی بین متفیرهای سیل 

شاخص‌های عامل تورم واریانس (VIF!)‏ و تولرانس CETOL)‏ به‌منظور برآورد همبستگی بین متغیرهای سیل 
استفاده گردید. مقادیر VIF‏ بزرگ‌تر از ۱۰ و مقادیر TOL‏ کمتر از ۰/۱ نشان‌دهنده هم خطی شدید است (منارد" 
É pla ۱‏ و همکاران ۲۰۰۵). 


TOL = 1-7 )۱( رابطه‎ 
VIF = 1, )۲( رابطه‎ 


که در آن RG‏ مقدار ضریب تعبین چندگانه بین متغیرهاست. مقدار ضریب تولرانس بین ۰ و ۱ متغیر است. هرچه 
مقدار آن بزرگ‌تر باشد (نزدیک به ۱)» میزان همپوشی با متغیرهای دیگر و در نتیجه هم خطی کمتر است. عامل تورم 
واریانس با تولرانس نسبت معکوس دارد. 

گام سوم: ایجاد مجموعه داده‌های آموزشی و اعتبارسنجی 

در این مرحله به‌منظور ایجاد مجموعه داده‌های آموزشی و اعتبارسنجی از روش نمونه‌برداری استفاده گردید. 
بدین ترتیب با استفاده از دستگاه GPS‏ و بازدیدهای میدانی و بهره‌گیری از امکانات سامانه Google Earth‏ و همچنین 
اطلاعات دبی و بارش شرکت مدیریت منابع آب ایران مناطق درگیر با حطر سیل گرفتگی (۱۹۰ نقطه) و مناطق بدون 
سیل (۱۹۰ نقطه) angi‏ شده و وارد پایگاه داده در محیط نرم‌افزار ArcGIS‏ شدند. سپس از این مجموعه داده‌ها؛ ۷۰ 
درصد برای آموزش (۱۳۳ نقطه برای سیل گرفتگی و ۱۳۳ مورد برای نقاط بدون سیل) و ۲۰ درصد OV)‏ نقطه برای 


سیل گرفتگی و ۵۷ مورد برای ble‏ بدون سیل) جهت اعتبارسنجی نتایج مدل‌ها مورد استفاده قرار گرفت. 


1 variance inflation factor 
2 tolerance 

3 Menard 

4 Hair 
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گام چهارم: مدل FURIA‏ 

FURIA‏ یک روش طبقه‌بندی فازی مبتنی بر قواعد است. این الگوریتم بر مبنای RIPPER‏ که یک الگوریتم 
یاد گیرنده قواعد یقینی است (کوهن '» ۱۹۹۵) توسعه یافته است. این روش توسط هان و هالرمیر" در سال ۲۰۰۹ 
معرفی شده است. FURIA‏ قوانین فازی را به جای قوانین مرسوم و مجموعه قوانین بی‌نظم را به جای فهرست 
قوانین یاد می‌دهد. این مدل به‌عنوان یک مدل طبقه‌بندی بسیار دقیق به نظر می‌رسد. که امکان تخمین مرزهای 
رشد و هرس کردن است. هرس کردن تا زمانی که میزان خطا بیشتر از ۵۰ درصد باشد ادامه می‌یابد (تراوینسکی " و 
همکاران ۲۰۱۱). به عبارت ساده‌تر قوانینی که تأثیری بر روی دقت مدل ندارد ضعیف شناخته شده و حذف 
می‌شوند. FURIA‏ بازه‌های یقینی به دست آمده از الگوریتم RIPPER‏ را به قوانین فازی که از توابع عضویت 
ذوزنقه‌ای استفاده می eS‏ تبدیل می کند (هان و Jla‏ مب ۰ در نتیجه قوانین فازی بر اساس PS a3‏ مقدار ۷ 


برای کلاس Yi‏ و درحه اطمینان CDi‏ ساخته می‌شوند: 
m‏ 
V = ( trate ((V) x CD; )۳( daly‏ 
i=1‏ 


گام پنجم: الگوریتم ژنتیک 

برای روش جستجوی استفاده شده توسط ارزیاب کننده FURIA‏ از الگوریتم ژنتیک استفاده شده است. 
الگوریتم ژنتیک یک روش فرا اکتشافی است که توسط هالند در سال ۱۹۷۵ ارائه شده است. GA‏ یک فرآیند مشابه 
برای انتخاب طبیعی به کار می‌گیرد و عمدتاً به‌عنوان یک روش جستجو برای یافتن راه‌حل‌های بهینه از طریق 
نظرسنجی از راه‌حل‌ه ای احتمالی استفاده می‌شود (کازوقلو" و همکاران» ۲۰۱۵). انتخاب تقاطع و جهش از 
عملگرهای پایه‌ای الگویتم ژنتیک برای تکثیر هستند (میقل CEN‏ این فرآیند تکثیر بهترین اعضاء که با عملگر 
هزینه ارزیابی می‌شونده به نسل بعدی منتقل می‌شوند. این فرآیند زمانی به پایان می‌رسد که آستانه‌های حاصی برای 


ارائه بهترین راه‌حل به دست می‌آیند (هاپت و Seale‏ ۲۰۰۶). 


1 Cohen 

2 Hiihn and Hiillermeier 
3 Trawinski 

4 Kavzoglu 

5 Mitchell 

6 Haupt and Haupt 


گام ششم: الگوریتم یادگیری کوانتیزاسیون برداری LVQ!)‏ 

بعد از انتخاب متغیرهای پیش‌بینی کننده» جهت بررسی رابطه بین خطر سیل با متغیرهای انتخاب شده از روش 
الگوریتم یادگیری کوانتیزاسیون برداری (1.۷60) استفاده گردید. 1.۷0۵ یکی از انواع شبکه‌های عصبی با الگوی 
ناد گیری نظارت‌شده است که در سال ۱۹۹٩‏ برای اولین بار توسط کوهن ارائه شده است. این روش یک رویکرد 
یادگیری الگو را به کار می‌برد که در واقع فاصله اقلیدسی را به‌عنوان یک قانون اساسی رقابت در نظر می گیرد. 

گام هفتم: روش ترکیبی LogitBoost‏ 

این روش را فریدمان و همکاران در سال ۲۰۰۰ معرفی کردند (فریدمان" و همکاران ۲۰۰۰). این الگوریتم 
طبقه‌بندی را با استفاده از یک برنامه رگرسیون به عنوان یادگیرنده پایه انجام می‌دهد و می‌تواند مشکلات چند کلاسه 
را Cs pte‏ کند (weka.classifiers.meta.LogitBoost)‏ در تحقیق حاضر درخت تصمیم گیری به‌عنوان طبقه‌بندی 
کننده پایه استفاده شده است (ایبا و لانگلی" ۱۹۹۲). در این پژوهش به‌منظور اجرای مدل ترکیبی از نرم‌افزار 
داده‌کاوی WEKA‏ استفاده شده است. نرم‌افزار WEKA‏ به عنوان یک نرم‌افزار یادگیری ماشین برای نخستین بار در 
کشور نیوزیلند و در دانشگاه Waikato‏ معرفی شده است. نام این نرم‌افزار از حروف اول عبارت Waikato‏ 
Environment for Knowledge Analysis‏ گرفته شده است. این نرم‌افزار مجموعه‌ای از به‌روزترین الگوریتم‌های 
یادگیری ماشین را جهت انجام داده‌کاوی ارائه می‌دهد. WEKA‏ شامل ابزارهایی برای آماده‌سازی داده‌ها؛ طبقه‌بندی؛ 
رگرسیون, خوشه‌بندی ایجاد قوانین انجمنی و تصویرسازی (کار با نمودارهای مختلف) است. 

گام هشتم: ارزیابی دقت نقشه حساسیت خطر وقوع سیل 

در پژوهش حاضر به‌منظور ارزیابی دقت نقشه حساسیت خطر وقوع سیل از منحنی مشخصه عملیاتی 
دریافت‌کننده یا منحنی مشخصه عملکرد سیستم (ROC)‏ و سطح زیر منحنی (AUC?)‏ استفاده شده است. در منحنی 
ROC‏ روی محور X‏ مقدار تشخیص‌پذیری یا ویژگی Specificity‏ (نسبت پیکسل‌های بدون سیل که به درستی 
به‌عنوان بدون سیل طبقه‌بندی شده‌اند) و روی محور ۷ مقدار حساسیت Sensitivity‏ (نسبت پیکسل‌های سیل که به 
درستی به‌عنوان سیل طبقه‌بندی شده‌اند) مشخص می‌باشد. مقادیر Specificity‏ و Sensitivity‏ با استفاده از روابط زیر 


محاسبه می‌شوند: 


1 Learning Vector Quantization 

2 Friedman 

3 Iba and Langley 

4 Receiver Operating Characteristics 
5 Area under the ROC curve 
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TN 
specificity = ————— (£) رابطه‎ 
P Y FP eai 
sensitivity = ۳ )۵( رابطه‎ 


که در آن TNI‏ تعداد مناطق بدون سیل که درست تشخیص داده شده‌اند؛ ۳۳2 تعداد تقاط سیل که اشتباه 
تشخیص داده شده‌اند؛ "1۳ تعداد نقاط سیل که درست تشخیص oala‏ شده‌اند و FNE‏ تعداد نقاط بدون سیل که اشتباه 
تشخیص داده شده‌اند. جهت نشان دادن عملکرد هر یک از مدل‌های ترکیبی از شاخص آماری صحت (Accuracy)‏ 
نیز استفاده شده است. این شاخص نسبت نقاط سیل گیر و بدون سیل که به درستی طبقه‌بندی شده‌اند را نشان 
می‌دهد. 


TP + TN 


(CV) رابطه‎ 
TP + TN + FP + FN 


Accuracy = 


منحنی ROC‏ در زمینه‌های مختلفی از جمله مهندسی» پزشکی و فیزیک به کار برده شده اما به طور عمده برای 
ارزیابی نقشه‌های پیش‌بینی مخاطرات طبیعی استفاده می‌شود (لی و همکاران» ۲۰۱۲؛ پرادهان و همکاران ٤۲۰۱؛‏ 
روی و ساها"؛ ۲۰۱۹). سطح زیر منحنی ROC‏ نشان‌دهنده توانایی یک مدل برای پیش‌بینی پیکسل‌های سیل و بدون 
سیل می‌باشد. مقدار AUC‏ صفر بیانگر این است که مدل غیرواقعی می‌باشد و هر چقدر مقدار AUC‏ بیشتر باشد 


عملکرد مدل در پیش‌بینی بهتر خواهد بود (تين بوی و همکاران. (YONA‏ 


۳- بحث و نتایج 


به‌منظور تهیه نقشه حساسیت خطر وقوع سیل در سطح حوضه آبریز الندچای ۱۳ پارامتر هیدروژئومورفیک مورد 
تجزیه‌وتحلیل قرار گرفت. 

لیتولوژی: خاک یا سنگ قابل‌نفوذ, شرایط نفوذ آب به داخل زمین را فراهم می‌کند و تخلیه آن را به داخل آبراهه 
اصلی به تأخیر می‌اندازد؛ ازاین‌رو رواناب سطحی کاهش می‌یابد. حوضه‌های با سنگ بستر یا خاک به نسبت غیرقابل 
نفو حجم بالایی از رواناب سطحی ایجاد می‌کنند (گرید ۲۰۰۹). نقشه لیتولوژی منطقه با استفاده از نقشه‌های 


1 True Negative 
2 False Positive 
3 True Positive 
4 False Negative 
5 Roy and Saha 
6 Garde 
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محیط نرم‌افزار ArcGIS‏ تهیه گردید. بررسی نقشه لیتولوژی OLS‏ می‌دهد که مناطق کم ارتفاع در قسمت شرقی 
منطقه را phy‏ سازندهای dle‏ مثل پادگانه‌های آبرفتی مربوط به دوره کواترنر و کنگلومرا تشکیل می‌دهد. عمده 
فعالیت‌های انسانی از قبیل شهرسازی و زمین‌های کشاورزی در این بخش از منطقه متمرکز شده‌اند. در مقابل 
ارتفاعات منطقه که بیشتر قسمت غرب و شمال غرب منطقه را دربرگرفته‌اند از سنگ‌های آتشفشانی بازالتی و در 
بخش‌هایی همراه با سنگ Sal‏ و شیل تشکیل شده‌اند. به دلیل شیب تقریباً زیاد این بخش‌ها و نفوذپذیری پایین 
سازندهاء آب‌های حاصل از بارندگی به‌صورت رواناب‌های متعدد شکل گرفته و باعث افزایش ناگهانی دبی رودخانه 


الندجای می‌شوند و در نتیجه در پایین‌دست حوضه» سیل گیری زمین‌های کشاورزی و مناطق مسکونی رابه همراه 


دارد. 
44°30'0"E 44°45'0"E‏ 
Z z‏ 
a =‏ 
Š ۳‏ 
a x‏ 
Loe]‏ ۹ 
آبرفت های بستر رودخانه ای 
پاد گانه های آبرفتی قدیم 

3 پاد گانه های آبرفتی جدید z‏ 
2 | سنگ آهک نومولیت دار S|‏ 
[A‏ سنگ آهک بلاژ یک بت 2 
2% ا ماسه سنگ همراه با کنگلومراء مارل و شیل 8 
۰ | گدازه های بالشی بازالتی خاکستری متمایل به سبز کنگلومرا 

سنگ مرمر آندزیت 

سنگ cle‏ آتفشانی بازالتی شیل و سنگ آهک 

ey‏ و سنگ آهک 
دیوریت و گنیس 
گنیس 
دخانه الندجا 

آمفیبولیت: ay‏ و دیوریت رو چای “NS,‏ 

گابرو و دیوریت 

سنگ های اولتراباز یک 18 12 6 3 0 

Ob >‏ های گدازه بازالتی همراه با سنگ آهک و شیل e‏ ات 

44°30'0"E 44°45'0"E 


شکل ۲- نقشه لیتولوژی منطقه مورد مطالعه 
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گروه‌های هیدرولوژیکی خاک: این پارامتر کیفیت SLE‏ را بر اساس حداقل میزان نفوذ آب نشان می‌دهند 
(دپارتمان کشاورزی ایالات متحده (VW‏ گروه‌های هیدرولوژیکی خاک به چهار گروه ۵ 8 ° و D‏ طبقه‌بندی 
می‌شوند. خاک‌های گروه A‏ دارای حداقل پتانسیل جهت تولید رواناب هستند؛ درحالی که خاک‌های گروه D‏ دارای 
حداکثر پتانسیل برای ایجاد رواناب می‌باشند (گیتلمان و همکاران؛ ۲۰۱۷). جهت ag‏ نقشه گروه‌های هیدرولوژیکی 
خاک eel plas‏ از نها زین شناسی> کاربری اراضی و همین داده‌های (Gling 8 She‏ 
هیدرلوژیک خاک برگرفته از وبسایت ناسا استفاده گردید. نقشه گروه‌های هیدرولوژیکی خاک منطقه در سه گروه B‏ 
C‏ و D‏ تهیه شده است. Bog S‏ زمین‌هایی را شامل می‌شود که از پادگانه‌ها و آبرفت‌های کواترنری پوشیده شده‌اند 
که به جهت نفوذپذیری نسبتاً حوب از پتانسیل کمی جهت تولید رواناب برخوردارند. گروه C‏ که بخش عمده‌ای از 
منطقه را شامل می‌شود اکثراً زمین‌های مرتفع و کوهستانی را در بر گرفته است. این مناطق به دلیل داشتن سازندهای 
نفوذناپذیر و همچنین رخنمون‌های سنگی از پتانسیل نسبتاً بالایی جهت تولید رواناب نسبت به گروه 8 برخوردار 
می‌باشند. گروه D‏ که پتانسیل بالایی را جهت تولید رواناب نشان می‌دهد. مناطق شهری و روستایی را شامل می‌شود. 
این مناطق به دلیل داشتن پوشش غیرقابل نفوذ از قبیل آسفالت و بتن شرایط مناسبی را جهت تولید رواناب فراهم 
می‌آورند. 
شاخص پوشش گیاهی: NDVI‏ شاخصی است که خصوصیات پوشش گیاهی یک منطقه را توصیف می کند. 
Shee peor‏ که بر رواناب سطحی و قابلیت نفوذپذیری یک منطقه تأثیر می‌گذارد (سویک و توپال" ۲۰۰۳؛ تهرانی» 
۳ مقدار آب و بار رسوب. پیک‌ه ای سیلاب و زمان وقوع آنه و سرعت انتقال پیک جریان به‌شدت توسط 
ماهیت و گسترش پوشش گیاهی تحت SE‏ قرار می‌گیرد (گرید» ۲۰۰۳). به‌منظور تهیه نقشه تراکم پوشش ALS‏ 
منطقه از شاخحص تفاضلی پوشش گیاهی نرمال شده استفاده شده است. 
کاربری اراضی: فراوانی وقوع سیل در یک منطقه می‌تواند به‌شدت تحت تأثیر الگوهای کاربری زمین و سیر 
تحول زمانی ol‏ قرار گیرد (بنیتو و همکاران" ۲۰۱۰؛ گارسیا رویز و همکاران؛ ۲۰۰۸). بیکرز و همکاران AVON)‏ 
در تحقیق خود نشان دادند که تغبیرات در کاربری اراضی می‌تواند احتمال وقوع سیل در یک منطقه را سرعت بخشد. 
جت تھ iat les aie‏ ری ماهوا اى ك ا سج shay by SOU‏ داد ال ۱۳۹۵ 
استفاده شده است. به همین منظور ابتدا پیش‌پردازش‌های لازم همچون تصحیح اتمسفری و رادیومتریکی روی 
تصوير اعمال شد. در مرحله بعد با استفاده از روش طبقه‌بندی نظارت‌شده و الگوریتم حداکثر احتمال نقشه کاربری 
ازاف درا R‏ ال ی کت 
Gittleman et al‏ 1 
Cevik and Topal‏ 2 


3 Benito et al 
4 Garcia-Ruiz et al 


سال دوازدهم مدل‌سازی حساسیت خطر وقوع سیل در حوضه آبریز ng‏ ۱۱ 


ارتفاع: ارتفاع یکی از مهم‌ترین متغیرهای تأثیرگذار بر روی سیل می‌باشد (تهرانی و همکاران» ۲۰۱6؛ بوی و 
همکاران (VENT‏ به‌طورکلی رابطه معکوسی بین خطر وقوع سیل و ارتفاع وجود دارد. فراوانی وقوع سیل با اف زایش 
ارتفاع کاهش bh ge‏ در نتیجه ارتفاعات پست‌تر از حساسیت بیشتری برای وقوع سیل برخوردار هستند (خسروی و 
همکاران ۲۰۱7). نقشه طبقات ارتفاعی منطقه با استفاده مدل رقومی ارتفاعی (DEM)‏ با قدرت تفکیک مکانی ۱۲/۵ 
متر در ۵ کلاس تهیه شده است. مقادیر ارتفاعی منطقه بین ۱۰۹۳ تا ۳۹۳۸ متر ارتفاع از سطح دریا متغیر است. با 
توجه به اینکه بیش از 7۰ درصد از کل مساحت منطقه را نواحی کوهستانی و ارتفاعات AL‏ ۲۰۰۰ متر تشکیل 
می‌دهد. سرعت شکل GS‏ رواناب در این بخش‌ها بالا بوده و در نتیجه احتمال he‏ سیل گیری مناطق پایین‌دست 


حوضه بعد از هر بارندگی افزایش می‌باید. 


شاخص پوشش ALS‏ گروه های هیدرولوژیکی خاک 
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شکل ۳- نقشه‌های گروه‌های هید رولوژیکی Se‏ شاخص پوشش گیاهی. ارتفاع و کاربری اراضی 


1 Bui et al 


۱۲ جغرافیا و مخاطرات محیطی شمارة اول 


شیب: بین شیب یک منطقه و سرعت جریان سطحی یک همبستگی قوی مثبت می‌توان یافت (فرناندز و لوتز؛ 
۰ داس» ۲۰۱۸). علاوه بر این شیب فرآیند نفوذ را کنترل می‌کند. رواناب سطحی با افزایش شیب به‌طور قابل 
توجهی افزایش می‌باید و در نتیجه مقدار نفوذ کاهش می‌یابد (داس و پردشی" ۲۰۱۸). 

جهت شیب: یکی دیگر از پارامترهای اصلی مورفومتریک که با رابطه معکوس در ارتباط با وقوع سیل در نظر 
گرفته می‌شود. جهت شیب است. دلیل آن این است که جهت شیب می تواند شرایط هیدرولوژیکی یک منطقه را 
تحت SE‏ قرار دهد (ارکان اوغلو و گوکچوغلو" ۲۰۰۲؛ کاواک و کوندو* ۲۰۰۸). جهت شیب به‌طور غیرمستقيم در 
وقوع سیل تأثیر دارد. به‌عنوان مثال مناطق شیب ‌دار یک منطقه که در قسمت سایه قرار دارند. رطوست خاک نسبتا 
بالایی داشته و در نتیجه با رواناب زیاد همراه هستند (توفیق لو اسلام و همکاران» ۲۰۲۱). نقشه جهت شیب منطقه در 
نرم‌افزار ArcGIS‏ و با استفاده از لایه DEM‏ تهیه شده است. 

فاصله از آبراهه: جریان‌های رودخانه‌ای مسیرهای اصلی برای دبی سیلاب هستند و مناطقی که Koay‏ 
رودخانه‌ها قرار دارند مستعد سیل‌گیری هستند (آپرمان و همکاران؟ ۰۹ فاصله از رودخان» یک عامل مهم 
ژئومورفیک است که جهت تهیه دقیق نقشه خطر وقوع سیل باید در نظر گرفته شود. با افزايش فاصله» شیب و ارتفاع 
بیشتر می‌شود. در نتیجه مناطق دور از JUS‏ رودخانه آسیب‌پذیری کمتری از وقوع سیل دارند (داس» ۲۰۱۸). 

تراکم آپراهه: زمانی که در یک حوضه آبریز بارندگی اتفاق می‌افند. تراکم زهکشی بر روی جریان‌های اوج تأثیر 
می‌گذارد (چوبین و همکاران» ۲۰۱۹). کومار و همکاران (۲۰۱۷) در پژوهش خود OLE‏ دادند که در مناطق با تراکم 
زهکشی VE‏ در مقایسه با مناطقی که از تراکم زهکشی پایینی برخوردار هستند. رواناب سطحی بالاتری تشکیل 


می‌شود. 


1 Fernandez and Lutz 

2 Das and Pardeshi 

3 Ercanoglu and Gokceoglu 
4 Kwak and Kondoh 

5 Opperman et al 


سال دوازدهم مدل‌سازی حساسیت خطر وقوع سیل در حوضه آبریز a‏ ۳ 
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شکل 6- نقشه‌های شیب. جهت شیب. فاصله از آبراهه و تراکم آبراهه 


بارش: بارندگی رابطه معنی‌داری با دبی رودخانه داشته و به‌طور مستقیم وقوع سیل را کتترل می‌کند. بارش‌های 
غیرمنتظره در مناطق نیمه‌حشک می‌توانند وضعیت سیلابی ایجاد کنند (داس» ۲۰۱۸). بنابراین می‌توان گفت که یک 
بارندگی شدید در بازه زمانی کوتاه می‌تواند Cok‏ وقوع سیل گردد. به‌منظور تهیه نقشه بارش منطفه از داده‌های 
ایستگاه‌های باران‌سنجی برای سال‌های دارای اطلاعات بارش و ایستگاه سینوپتیک خوی برای یک دوره ۳۰ ساله 
(۱۳۹۹-۱۳۷۰) استفاده شده است. 

شاخص رطوبت تویوگرافیک: این پارامتر در سال ۱۹۷۹ توسط بیون و کرکبی ارائه شده است و نشان‌دهنده 
Ol‏ فضایی رطوبت در حوضه آبریز است (رحمتی و همکاران» ۲۰۱7). به عبارت دیگ شاخص رطوبت 
تویوگرافیک مقدار انباشت آب در هر پیکسل از حوضه آبربز را نشان می‌دهد (گوکچوغلو و همکاران ۲۰۰۵). 
IS) be‏ مقادیر بالای شاخص رطوبت توپوگرافیک و وقوع سیل همبستگی قوی با SAS‏ دارند (شیت و 
همکاران ؛ ۲۰۲۰). برعکس, مناطقی که polia‏ کمتری از این شاخص را نشان می‌دهند از آسیب‌پذیری کمتری در 


برابر سیل برخوردار هستند. 


1 Shit et al 


۱ جغرافیا و مخاطرات محیطی شمارة اول 


شاخص قدرت آبراهه: Gb‏ نظر نایتون" )1994( شاخحص قدرت آبراهه برای بسیاری از فرآیندها در محیط‌های 
رودخانه‌ای از اهمیت بسیار بالایی برخوردار است. این شاخص را می‌توان به‌عنوان یک عملکرد قابل توجه در 
فرسایش کانال رودخانه و حمل رسوب توصیف کرد Sp)‏ و همکاران"؛ ۲۰۰۹؛ هانگ و همکاران ۲۰۱۸). 

شاخص حمل رسوب: این شاخص نشان‌دهنده قدرت جریان آب از نظر فرسایش است که بر روی شرایط 
میدرولوژیکی اثر می‌گذارد (تهرانی و همکاران» ۲۰۱۵). هانگ و همکاران (۲۰۱۸) در پژوهش خود به این نتیجه 
رسیده‌اند مناطقی که مقادیر STI‏ کمتری دارند اکثراً در زمین‌های مسطح قرار دارند و بنابراین پتانسیل بالایی برای 


شاخص رطوبت توپ و گرافیکك 
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شکل ۵- نقشه‌های بارش. شاخص رطوبت تویوگرافیک. شاخص قدرت آبراهه و شاخص حمل رسوب‎ 


1 Knighton 
2 Barker et al 


سال دوازدهم مدل‌سازی حساسیت خطر وقوع سیل در حوضه آبریز .... ۱۵ 


تحلیل عوامل هم خطی سیلاب 

در این بخش از تحقیق تجزیه‌وتحلیل هم خحطی چندگانه به‌منظور بررسی همبستگی داخلی متغیرها انجام شده 
است که نتایج آن نشان می‌دهد تقریباً بین هیچ یک از پارامترهای مرتبط با سیل هم خطی چندگانه جدی وجود ندارد 
(جدول ۱). پارامتر گروه‌های هیدرولوژیکی خاک کم‌ترین مقدار تولرانس را داشته است (۰/۱۳۹)» که بیشتر از مقدار 
بحرانی نظری (۰/۱۰۰) می‌باشد. پارامترهای شاخص حمل رسوب و شیب در رتبه بعدی قرار دارند. همچنین برای 
هر یک از پارامترها مقدار عامل تورم واریانس نیز محاسبه گردید که نتایج آن نشان داد پارامتر گروه‌های 
هیدرولوژیکی خاک با وجود داشتن بالاترین مقدار یعنی ۷/۱۹ از حد آستانه نظری برای هم حطی چندگانه خیلی 
فاصله دارد (بیشتر از ۱۰). سایر پارامترها مقادیر پایین‌تری را نسبت به گروه‌های هیدرولوژیکی خاک نشان می‌دهند. 
بنابراین همه پارامترها در ارتباط با حساسیت خطر وقوع سیل در نظر گرفته شده و در گام بعدی مورد تجزیه‌وتحلیل 


قرار گرفته‌اند. 


پارامترهای مرتبط با سیل Variance IF (VIF) Tolerance (TOL)‏ 
لیتولوژی ۷۰۱ VEY‏ 
گروه‌های هیدرولوژیکی خاک Ya‏ ۷/۹ 
پوشش گیاهی “YY‏ 1/1۹ 
کاربری اراضی YY YY)‏ 
ارتفاع v/s) YY‏ 
شيت Yt‏ 11۹۰ 
جهت شیب 7۱ ۱/۳۹ 
فاصله از آبراهه ۳۳۲ ve)‏ 
تراکم آبراهه abi ‘AAs‏ 
بارش ۰/۳۹ ۳/۳ 
شاخحص رطوبت توپوگرافیک ۰۹۹ ۱/۳۸ 
شاخحص قدرت آبراهه ۳۹ ۳/۳ 
شاخص حمل رسوب ۰/4 5/۹ 


۱۳1 جغرافیا و مخاطرات محیطی شمارة اول 


اجرای الگوریتم FURIA-GA‏ 


جهت اجرای الگوریتم FURIA-GA‏ خصوصیات زیر که پس از انجام آزمون و خطای گسترده به دست 
آمده‌اند. تعیین شدند. برای GA‏ احتمال هم‌گذری' بر روی ۲ احتمال her‏ روی ۵ اندازه جمعیت" 
۰ و تعداد نسل‌ها بر روی ۰ تنظیم شد. برای ارزیاب کننده FURIA‏ یک تکنیک اعتبارسنجی متقاطع ۱۰ 
برابری" با product T-Norm‏ به‌عنوان عملگر اجتماع فازی به‌منظور تر کیب فالخ شین اموزشن داده شد. نتایج 
نشان داد که طبقه‌بندی FURIA-GA‏ با مقدار AVEO‏ / دقت خیلی خوبی را داشته است. 


تعیین اهمیت متغیرها 
از روش LVO‏ جهت تعیین اهمیت پارامترهای مرتبط با سیل که در مرحله قبل با استفاده از الگوریتم FURIA-‏ 
GA‏ انتخاب let‏ در نرم‌افزار WEKA‏ استفاده گردید. مقادیر LVQ‏ بین EA‏ ۳ ۰" متغیر بوده است که 


بالاترین مقدار مربوط به شیب و کمترین مقدار نیز مربوط به شاحص رطوبت توپوگرافیک می‌باشد (شکل AV‏ 


شاخص رطوبت توپوگرافیک 
شاخص حمل رسوب 

بارش 

تراکم آبراهه 

لیتولوژی 

فاصله از آپراهه 

شاخص قدرت آپراهه 
کاربری اراضی 

پوشش گیاهی 

ارتفاع 

گروه های هیدرولوژیکی خاک 


شکل -I‏ اهمیت پارامترهای مرتبط با سیل بر اساس روش LVQ‏ 


1 crossover probability 

2 mutation probability 

3 population size 

4 number of generations 

5 fold cross-validation technique 


سال دوازدهم مدل‌سازی حساسیت خطر وقوع سیل در حوضه آبریز a‏ ۷ 


احرای روش ترکیبی LogitBoost‏ 

به‌منظور اجرای الگوریتم LogitBoost‏ در نرم‌افزار WEKA‏ تنظیمات زیر استفاده شد: اندازه بسته ؛ ٩۱۰۰‏ تعداد 
TLS‏ ۱۲؛ تعداد بذر یا Tails‏ ۱. همچنین decision stump tree‏ به‌عنوان طبقه‌بندی کننده پایه" انتخاب گردید. با 
توجه به اینکه فرمت اصلی داده‌های WEKA‏ به‌صورت ARFF‏ می‌باشد. به‌منظور انتقال لایه‌های اطلاعاتی به محیط 
نرم‌افزار WEKA‏ و اجرای هر یک از مراحل فوق» ابتدا این لایه‌ها به فرمت ASCH‏ تبدیل شده و سپس با استفاده از 
نرم‌افز ار SPSS‏ در محیط WEKA‏ بارگذاری شدند. 


angi‏ نقشه حساسیت خطر وقوع سیل 
در این بخش از تحقیق به‌منظور تهیه نقشه حساسیت خطر وقوع سیل» نتایج حاصل از اجرای مراحل ذکر شده 
در روش تحقیق در نرم‌افزار WEKA‏ جهت تبدیل به نقشه مصور به نرم‌افزار ArcGIS‏ منتقل شدند. نقشه حساسیت 
خطر وقوع سیل در سطح حوضه آبریز الندچای با استفاده از روش Natural Break‏ در ۵ کلاس تحت عنوان 
حساسیت خیلی کم. کم متوسط زیاد و خیلی زیاد تهیه گردید (شکل (V‏ نقشه نهایی به‌صورت پیکسل پایه بوده و 
حساسیت هر پیکسل را نسبت به خطر وقوع سیل نشان می‌دهد. 
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شکل ۷- نقشه نهایی حساسیت خطر وقوع سیل در حوضه ابریز الندجای‎ 


1 Batchsize 

2 number of iterations 
3 seed 

4 base classifier 


۱/۸ جغرافیا و مخاطرات محیطی شمارة اول 


بررسی نقشه نهایی نشان می‌دهد مناطقی که از نظر خطر وقوع سیل حساسیت بالایی را نشان می‌دهند عمدت 
قسمت‌های شرقی و جنوب شرقی حوضه را شامل می‌شوند. این قسمت‌ها عمدتاً در مناطق پایین‌دست حوضه قرار 
گرفته‌اند که هم ارتفاع حیلی کم و هم شیب کمی نسبت به ناهمواری‌های اطراف خود دارند. بنابراین رواناب‌هایی که 
پس از هر بارندگی یا با ذوب برف‌ها شکل می‌گیرند در این مناطق متمرکز می‌شوند که باعث سیل گرفتگی مناطق 
شهری و روستایی و زمین‌های کشاورزی و باغات می‌شود. مهم‌ترین مرکز جمعیتی حوضه یعنی شهر خوی نیز در 
این قسمت جای گرفته است که آسیب‌پذیری بالای تأسیسات شهری را هنگام وقوع سیلاب‌های شدید و ناگهانی 
نشان می‌دهد. به‌طورکلی روستاهایی که در مسیر رودخانه الندچای شکل گرفته‌اند از آسیب‌پذیری بالایی هنگام وقوع 
سیل برخوردار هستند. پهنه‌هایی که در کلاس‌های زیاد و خیلی زیاد از نظر حساسیت خطر وقوع سیل قرار دارند 
بیش از ۲۰ درصد از کل مساحت منطقه را شامل می‌شوند. بیشترین مساحت نیز مربوط به کلاس کم می‌باشد که با 
مقدار ۵۷۳/۱۲ کیلومترمربع نزدیک به ۵۰ درصد از مساحت کل منطقه را پوشش داده است. مساحت و درصد 


مساحت هر یک از کلاس‌های خطر وقوع سیل در جدول ۲ ارائه شده است. 


جدول ۲- مساحت و درصد مساحت هر یک از کلاس‌های خطر وقوع سیل 


کلاس مساحت (کیلومترمربع) درصد مساحت 
خیلی کم \v/W Ye YVA‏ 
کم ovy Y‏ £4/40 
متوسط 3۸1 3/۱ 
زیاد viv" Wy»‏ 
خیلی VEAL Lj‏ ۱۳/۳ 


ارزیابی عملکرد و دقت مدل 

هدف نهایی در agi‏ نقشه‌های حساسیت خطر وقوع See‏ پیش‌بینی مناطقی هست که احتمال دارد در آینده با 
وقوع بارش‌های شدید تحت تأثیر سیل قرار بگیرند. بنابراین بدون توجه به میزان دقت نقشه‌های نهایی تهیه شده با 
مدل‌ها و الگوریتم‌های متخلف استناد به این نقشه‌ها از اعتبار علمی لازم برخوردار نخواهد بود. لذا در این بخش از 
تحقیق به‌منظور ارزیابی عملکرد و دقت مدل مورد استفاده در تهیه نقشه نهایی حساسیت خطر وقوع سیل در سطح 
حوضه al‏ یز الندجای از ۳ شاخص آمار ی Accuracy Specificity Sensitivity‏ منحنی ROC‏ و همچنین سطح 
زیر منحنی در نرم‌افزار SPSS‏ استفاده شده است. نتایج بررسی ۳ شاخص آماری نشان می‌دهد که مدل FURIA-‏ 


GA LogitBoost‏ با ضریب ۸۷۰۹ از دقت خوبی بر اساس داده‌های آموزشی برخوردار می‌باشد (جدول (Y‏ منحنی 
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ROC‏ نیز برای دو گروه داده‌های آموزشی و اعتبارسنجی ترسیم شد. نتایج نشان داد که در مورد داده‌های آموزشی؛ 


مدل FURIA-GA-LogitBoost‏ با سطح زیر منحنی برابر با ۰/۸۲۱ از دقت خوبی برخوردار می‌باشد. 


جدول ۳- نتایج عملکرد مدل ترکیبی بر اساس داده‌های آموزشی و اعتبارسنجی 


داده‌های اعتبارسنجی داده‌های آموزشی شاخص آماری 
True positive ۰ oY‏ 
True negative 1٤ Os‏ 
False positive \A ۵‏ 
False negative ۹ ۷‏ 
Sensitivity (%) Ao/AY MWY‏ 
Specificity (%) AV 1‏ 
Accuracy (%) ۸7۰۹ AVE‏ 
O N E‏ 1 
0.8 
0.6 2 
= 
5 
c‏ 
v‏ 
N‏ 
0.4 
Training (AUC= 0.861)‏ —— 
0.2 
validation (AUC= 0.895)‏ -== 
Reference Line‏ — 
0.0 
1.0 0.8 0.6 0.4 0.2 0.0 
Specificity‏ - 1 


JSS‏ ۸- منحنی ROC‏ بر اساس داده‌های آموزشی و اعتبارسنجی 


۶- جمع‌بندی 

در تحقیق حاضر تلاش شد با استفاده از یک مدل ترکیبی نوین به مدل‌سازی حساسیت خطر وقوع سیل در 
حوضه آبریز الندچای پرداخته شود؛ بنابراین ابتدا ۱۳ پارامتر مؤثر جهت تهیه نقشه حساسیت خطر وقوع سیل انتضاب 
گردید. بررسی اهمیت پارامترها با استفاده از الگوریتم یادگیری کوانتیزاسیون برداری نشان داد که پارامتر شیب با 
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وقوع سیل در سطح حوضه آبریز الندچای در ۵ کلاس از خطر خیلی کم تا خیلی زیاد تهیه گردید که بررسی 
مساحت هر یک از کلاس‌ها نشان داد بیش از ۲۰ درصد از کل مساحت منطقه در پهنه‌های با حساسیت زیاد و خیلی 
زیاد قرار دارند. این پهنه‌ها که عمدتاً در پایین‌دست حوضه قرار دارند محل شکل‌گیری اجتماعات انسانی و 
زمین‌های کشاورزی هستند که با طغیانی شدن رودخانه الندچای خسارات زیادی به آن‌ها وارد می‌شود. در گام نهایی 
تحقیق میزان دقت و عملکرد مدل مورد استفاده در تحقیق با استفاده از منحنی ROC‏ مورد ارزیابی قرار گرفت که 
نتیجه Sly of‏ عملکرد خوب مدل FURIA-GA-LogitBoost‏ در تهیه نقشه حساسیت خطر وقوع سیل است. 
یافته‌های پژوهش حاضر با نتایج محققینی مثل تین بوی و همکاران (۲۰۱۹) در استان لائو سای واقع در کشور ویتنام 
مطابقت دارد. محققین مذکور در تحقیق خود به عملکرد خوب مدل‌های ترکیسی نوین جهت مدل‌سازی و تهیه 
نقشه‌های حساسیت خطر وقوع سیل اشاره می‌کنند که نتایج به‌دست‌آمده از پژوهش حاضر نیز بیانگر تأیید نتایج این 


قدردانی 
تحقیق حاضر برگرفته از نتایج رساله دکتری در دانشگاه تبریز است و از طرف صندوق حمایت از پژوهشگران و 
فناوران کشور (Iran National Science Foundation: INSF)‏ با کد ۹۸۰۲۷۱۵۲ حمایت مالی شده است. که بدین 


وسیله از آن‌ها قدردانی می‌شود. 
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